Thermodynamical properties of a trapped interacting Bose gas 



(N 

o 

(N 



Shi-Jie Yang*, Yuechan Liu, and Shiping Feng 
Department of Physics, Beijing Normal University, Beijing 100875, China 

The thermodynamical properties of interacting Bose atoms in a harmonic potential are studied 
within the mean-field approximation. For weak interactions, the quantum statistics is equivalent to 
an ideal gas in an effective mean-field potential. The eigenvalue of the Gross-Pitaevskii equation 
is identified as the chemical potential of the ideal gas. The condensation temperature and density 
profile of atoms are calculated. It is found that the critical temperature T c decreases as the inter- 
actions increase. Below the critical point, the condensation fraction exhibits a universal relation of 
No/N = 1 — (T/T c ) 7 , with the index 7 w 2.3 independent of the interaction strength, the chemical 
potential, as well as the frequency of the confining potential. 
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I. INTRODUCTION 

Since the first observations of the Bose-Einstein con- 
densation (BEC) in dilute alkali-metal atom gases [l|-[3j], 
experimental developments have posed many new tests 
for many-body theory. Numerous theoretical approaches 
have been implemented in order to obtain accurate re- 
sults for both the ground-state and non-equilibrium prop- 
erties of the trapped boson systems 043] ■ Although 
there is a vast number of studies in literature, the ef- 
fects of interactions on the transition temperature of 
BEC still remains an open question. The current con- 
sensus among physicists is that the interaction-caused 



shift in condensation temperature AT C = T c 



(«) 



T, 



(0) 



where tJ ^ and Tc are respectively the critical tem- 
perature for ideal and interacting gases, satisfies a lin- 
ear relation AT C /Tc ~ cagn 1 / 3 with respect to the s- 
wave scattering length a s . Here n is the particle den- 
sity. However, the value and even the sign of the coef- 
ficient c is controversial [8|. Some authors claimed that 
repulsive interactions reduce the critical temperature [9- 
Il2| . while others indicated an increase of the critical 
temperature [l3l - l2l1 | . In a recent publication p2|. the au- 
thors obtained a negative constant c « —2.33 by using 
the model of spatial permutations. This contradiction 
is somehow reconciled by a recent experiment [23}, where 
the researchers demonstrated a negative shift of the criti- 
cal temperature for weak interactions and a positive shift 
for sufficiently strong interactions. 

For an ideal Bose gas in a harmonic trap V(r) = 
^mcj 2 r 2 , the BEC temperature T c is determined by 

k B T c = hu{N/C(3)) 1/3 , with C(3) = 1.202(1 0, |H|. The 
fraction of condensation at a temperature T below T c 
is N /N = 1 - (T/T c ) 3 , where N is the number of 
Bose particles in the condensed state. In this paper, 
we explore the phase transition properties of interacting 
Bose atoms confined in a spherically symmetric harmonic 



potential within the mean- field approximation (MFA). 
We show that in the weak interaction limit the parti- 
cle distribution is equivalent to the ideal gas in an ef- 
fective mean-field potential. The eigenvalue of the gen- 
eralized Gross-Pitaevskii (GP) equation is identified as 
the chemical potential of the ideal gas. For the con- 
densed atoms we obtain the density profile within the 
Thomas-Fermi approximation. The condensation tem- 
perature T c decreases with increase of interactions as 
T c ~ etj^ with j3 sa —0.75, in contrast to the pre- 
vious mean-field calculations which gave a linear rela- 
tion AT C /T C (0) w -2.47a s n 1 / 3 [H|23. Our results qual- 
itatively coincide with the recent measurement (23| . In 
particular, the T c is independent of confining frequency 
which implies the details of the external potential is 
unimportant to the BEC. Below the critical tempera- 
ture, we find a universal relation between the condensed 
fraction and the temperature Nq/N = 1 — (T/T c ) 7 , with 
an index 7 w 2.3 being independent of the interaction 
strength, the chemical potential, as well as the frequency 
of the confining potential. 

The paper is organized as following: In Sec. II we briefly 
introduce the mean-field theory and the related results. 
Section III displays our main numerical results. A brief 
summary is included in Sec. IV. 



II. MEAN-FIELD THEORY 

The Hamiltonian for an interacting Bose gas in a ex- 
ternal trap V(r) is 



V}Vj- 



|h 4 }, 



(1) 
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where the coupling constant g — Ai:h 2 a s /m is written 
in terms of the scattering length a s and atom mass m. 
Following the standard approach, the Bose field operator 
is decomposed as tjj — $ + <f>, where $ is a c-number for 
the condensate, and </>(r) is an operator representing its 
fluctuations which annihilates a thermal atom at r. 
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Excluding the possibility of aggregate motion and vor- 
tices, the stationary state of $(r) satisfies the generalized 
GP equation[li-E3 

- |- V 2 $ + [V(r) - M + 9(n (r) + 2n T (r))]$ = 0, (2) 
2m 

where fi is the eigenvalue of the GP equation. The equa- 
tion of motion for (</)(r, t)) is 

ih ^m = -£v 2 (0) + ^(r)-M+2.9n(r)](0)+. 9 n o (r)(^). 

(3) 

The condensate density and the number of conden- 
sate atoms are defined by n (r) = | <&(!") | 2 and N = 
J G?rno(r), respectively. The noncondensate (thermal) 

density is nx — (4>'4>) : where (• • • ) indicates a thermal 
average in the grand-canonical ensemble. The general- 
ized nonlinear GP equation which includes interac- 
tions between the condensate and the thermal atoms, 
is solved for a static condensate, while the quasiparti- 
cle excitations of the system is described by the coupled 
Bogoliubov equations (O, which yields the quasiparticle 
energies and amplitudes. These in turn determine the 
number of noncondensed atoms Nt as well as various co- 
herence terms thermodynamic averages over two or more 
Bose field operators. 



As emphasized in Ref. 28|, the coherence terms = 
- $ 2 yield an excitation spectrum that is not gap- 
less. The Popov approximation, which neglects these 
terms, reduces the third and fourth-order terms in <j) and 
(ft, respectively, to the first and second order so that 
the Hamiltonian can be diagonalized 29] . Although it 
has been successful in describing the properties of the 
trapped Bose gases, it is not well-grounded theoretically, 
and fails to yield accurate predictions for the low-lying 
excitations at high temperatures. 

The elementary excitations in Eq.Q can be obtained 
by employing the semiclassical WKB approximation [l2j . 
The coupled equations for (<f>) and its complex conjugate 
(ify), can be solved explicitly, and the semiclassical exci- 
tation spectrum is obtained as 



e(p, r) = W + V(r) - /i + 2 5 n(r)] 2 - g 2 n 2 (r). (4) 

For a homogeneous system at low temperatures, tit can 
be neglected, and the above excitations coincide with the 
usual Bogoliubov spectrum. 

The quasiparticles with energies e(p, r) are distributed 
according to the Bose distribution function, 



/(p,r) = {exp[/3e(p,r)]-l}- 



(5) 



where /3 = l/fe^T. The particle distribution function can 
be obtained from the Bogoliubov canonical transforma- 
tions and is given by 



de , 



F(p,r) = -( — ) no /(p,r). 



(6) 



It is notable that the identification of the chemical po- 
tential ft with the eigenvalue /i of the GP equation is 
questionable in general. In the grand-canonical ensem- 
ble, the chemical potential is defined as ft = dE/dN, cor- 
responding to the energy cost E of adding a particle to 
the entire system, not only to the condensate. However, 
in the framework of the present theory, the excitation en- 
ergies (U) can be approximately reduced to the first order 
of coupling strength g as 



e(P,r) 



P 

2m 



+ V(r) - n + 2gn(r). 



(7) 



Therefore, (de/dfi) n ~ — 1, yielding the thermal density 
distribution 



n T (r) 



1 



dp 



exp[/3(p 2 /2m + V eff (r) - fi)] - 1 (2^)3 



(8) 



This equation is exactly equivalent to the statistics of an 
ideal Bose gas confined in an effective mean-field poten- 
tial, 



V eff = V{v) + 2gn{v). 



(9) 



Consequently, the eigenvalue fi of the GP equation Q is 
justified as the chemical potential of the ideal gas in the 
effective potential. 

The result can be understood as following: Suppose 
an interacting Bose gas is trapped in an external poten- 
tial V(r) which reaches an equilibrium distribution n(r). 
Adding an extra noninteracting Bose atom increases an 
energy Ae(r) = V(r) + 2gn(r). As the new comer is 
assimilated into the interacting gas, the density distribu- 
tion n(r) adjust slightly in order to transfer the increased 
energy into the effective potential. The factor 2 in the 
effective potential ([S]) comes from the Bose symmetry. 



III. NUMERICAL RESULTS 

At finite temperature T below T c , interaction effects 
involving the thermal atoms should be taken into account 
in addition to those of the condensate, especially when 
T goes close to T c . Integrating over the momentum p in 
Eq.©, we obtain in the harmonic potential, 



tit 



e -k/3(muj 2 r 2 /2+2gn(r)-fj.) 



k=l 



(10) 



where At = (2nh 2 /raksT) 1 / 2 is the thermal de Broglie 
wavelength. 

From the generalized GP equation @, the Thomas- 
Fermi approximation 



(j, = V(r)+g[n (r) + 2n T (r)], 



(11) 



will be used in the following discussions. For a spherically 
symmetric potential, there is a boundary radius r%. For 
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^ g °o ' 2 '4 i i — io FIG. 2: (Color online) The eondensed atoms versus the tem- 

f r perature for (a) various interaction strength at u) = 1 and 

fi = 6 and for (b) various chemical potentials at cj = 1 and 
3 = 0.5. 



FIG. 1: (Color online) Density distributions of the condensed 
(real curves) and noncondensed (dashed and dotted curves) 
atoms below the critical temperature, (a) Temperature de- 
pendence at fj, = 10, g = 0.1, and ui — 2. (b) Interaction 
strength dependence at n — 10, T = 3, and uj — 2. (c) Chem- 
ical potential dependence at g = 0.1, T = 3, and cj = 2. 
(d) Potential frequency dependence at fi = 10, T = 3, and 
3 = 0.1. 



r > ri, the condensate density no(r) = 0. The thermal 
density is determined by, 



n T (r) 



fe=i 



-fc/3[mw 2 r 2 /2+2gn T (r)-/i] 



(A T Vfc) E 



(12) 



Whereas for r < ri, %(r) 7^ 0, 

00 e k/3[mu> 2 r 2 /2+2gn T {r)~fj,] 



fc=l 



(A T Vfc) E 



(13) 



In Fig.l, we plot the various cases of the density pro- 
files of condensed and noncondensed atoms in the spher- 
ically symmetric harmonic potential. Hereafter the units 
of h = m = ks = 1 are used. For a better understand- 
ing of the dependence of the density on temperature, in 
Fig. 2 we plot the total condensed atoms as a function of 
temperature. Figure 2(a) shows that the number of con- 
densed atoms decreases with increase of interactions at 
a fixed chemical potential and a given temperature. Fig- 
ure 2(b) is the number of condensed atoms for different 
chemical potentials at a given interaction. It coincides 
with the measurements carried out by Ensher et al[30|. 
In particular, as expected, it increases with the chemical 
potential. 

We note that in our calculations that in our parameter 
ranges, the critical temperature is almost independent of 
the characteristic frequency w of the confining potential. 
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FIG. 3: (Color online) The condensed atoms versus the tem- 
perature for various frequencies of the confining potential at 
g = 0.1 and /i = 6. 



Figure 3 reveals that all curves meet at one temperature 
point where the condensed atoms tend to zero. We con- 
sider that this result may be relevant to the mean-field 
theory when applying to a finite system, as well as semi- 
classical approximation in the calculation of the state 
density. Nevertheless, it implies that the condensation 
temperature is insensitive to the details of the confining 
trap. 

At the critical temperature T c , the condensed atoms 
is nearly zero. We obtain the critical temperature from 
Fig. 2 for various interaction strength and chemical poten- 
tials, respectively, as shown in Fig. 4. Figure 4(a) shows 
that the critical temperature T c decreases with increase 
of the interaction strength g. The shift of critical tem- 
perature becomes saturated as the interactions increases. 
The inset is a log-log plot of the critical temperature ver- 
sus the interaction strength. It reveals that the critical 
temperature decreases with the interaction strength or s- 
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FIG. 4: (a) The critical temperature versus the interaction 
strength for ui — 1 and fi — 6. Inset: The Log-Log plot gives 
a nearly linear relation with a slope j3 ~ —0.75. (b) The 
critical temperature versus the chemical potential for to = 1 
and g — 0.5. 
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FIG. 5: (Color online) The log-log plot of the condensate 
fraction versus temperature for different interaction strength 
(a) and for different chemical potentials (b). 



wave scattering length a s as T c ~ af with j3 m —0.75, in 
contrast to previous MFA calculations which give a linear 
relation AT C /T C (0) w -2.47a s n 1 / 3 [H, H3). Our results is 
qualitatively consistent with the experimental measure- 
ments by Smith et al 23] . Although the mean-field theory 
has weakness in dealing with the thermodynamical prop- 



erties near the critical point, as discussed in Sec. II, it is 
remarkable that the results agree with the experiment 
quite well. Figure 4(b) shows that the critical tempera- 
ture almost linearly increases with the chemical potential 

To examine the universal behavior below the critical 
temperature, we study the temperature dependence of 
the condensate fraction Nq/N which takes the form, 



Figure 5 is a log-log plot of the condensate fraction ver- 
sus the reduced temperature. It shows that for vari- 
ous interaction strength (Fig. 5(a)) and chemical poten- 
tials (Fig. 5(b)), respectively, the lines are parallel to each 
other, giving rise to a common slope 7 ~ 2.3. This in- 
dex 7 is universal which is independent of the interaction 
strength, the chemical potential, as well as the frequency 
of the confining potential. In comparison, the ideal gas 
in a confining potential has an index 70 = 3, which is 
also independent of the characteristic frequency of the 
confining potential. 



IV. SUMMARY 

We have studied the thermodynamical properties of 
the interacting Bose gas confined in an external poten- 
tial. For weak interactions, the system can be equivalent 
to an ideal gas in an effective mean-field potential where 
the eigenvalue of the generalized GP equation plays the 
role of the chemical potential. The critical temperature 
decreases nonlinearly with increase of interactions. Our 
results agree quite well with the recent experimental mea- 
surements. Furthermore, we found a universal relation of 
the condensate fraction with respect to the temperature 
below the critical temperature, with the index being in- 
dependent of the interaction strength, the chemical po- 
tential, as well as the frequency of the confining poten- 
tial. We expect verifications of the universal relation by 
experiments. 
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